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ABSTRACT 


This  report  derives  from  elementary  principles  the  general  equations  of 
motion  for  a  missile  utilizing  a  fixed-plane  coordinate  system,  i.e.,  a  co¬ 
ordinate  system  with  one  axis  constrained  to  lie  in  a  given  plane. 

Included  in  the  derivation  are  explicit  expressions  for  introducing  v»:nd 
and  an  alternate  set  of  equations  to  be  used  when  singularity  conditions 
are  approached.  Means  are  provided  for  automatically  converting  to  the 
alternate  set  of  equations  so  that  uninterrupted  trajectory  simulation  can 
proceed  under  all  conditions.  A  complete  discussion  of  initial  conditions 
is  included. 

The  general  equations  can  be  used  for  flat  or  spherical,  rotating,  or  non¬ 
rotating  earth  cases. 


INTRODUCTION 


This  report  is  part  of  a  continuing  program  to  give  Picatinny  Arsenal 
a  complete  capability  in  the  flight  simulation  of  all  types  of  projectiles 
and  missiles,  whether  ballistic  or  rocket-boosted,  guided  or  unguided. 

Contained  herein  is  a  rederivation,  from  elementary  principles,  and  an 
elaboration  of  the  fixed  plane  coordinate  system  described  in  Reference  1, 
"Trajectory  Equations  for  a  Six-Degree-of-Freedom  Missile."  In  Refer¬ 
ence  1,  a  list  of  only  the  basic  equations  unique  to  the  fixed  plane  co- 
ordina  e  system  is  presented.  This  report  completes  this  list  and  includes, 
for  the  first  time,  a  derivation  for  initial  conditions,  as  well  as  an  alter¬ 
nate  set  of  equations  to  be  used  when  singularity  conditions  are  approached. 
Means  are  provided  for  converting  to  this  set  during  flight  simulation. 
Finally,  explicit  equations  for  introducing  wind  into  the  equations  of 
motion  are  derived. 

This  alternate  set  of  fixed-plane  trajectory  equations  was  required  to 
supplement  (and,  in  many  cases,  replace)  the  existing  equations  utilizing 
the  missile-fixed  coordinate  system  as  derived  in  Reference  1.  To  be 
explicit,  the  previous  equations  produced  satisfactory  trajectories  for  low- 
spin  projectiles,  but  were  inadequate  for  the  simulation  of  spin-stabilized 
shells.  The  most  obvious  differences  were  in  the  large  deflections  accom¬ 
panying  most  trajectories  for  high-spin  projectiles,  these  being  three  or 
four  times  as  large  as  had  been  anticipated.  Discussion  of  this  matter 
with  personnel  from  the  Naval  Weapons  Laboratory,  Dahlgren,  Virginia, 
indicated  the  existence  of  a  narrow  band  of  permissible  integration  step 
sizes  (incremental  time  steps)  whereby  both  truncation  and  round-off 
errors  are  of  acceptable  magnitudes.  It  is  possible  that  when  forces  and 
moments  are  referred  to  a  missile-fixed  coordinate  system  (as  was  the  case 
in  Reference  1),  and  high-spin  rates  arc  to  be  accounted  for,  this  band  of 
acceptable  time  increments  becomes  even  more  narrow  or  perhaps  non¬ 
existent.  Although  this  is  not  known  with  certainty,  it  provides  the  im¬ 
petus  for  studying  the  fixed-plane  coordinate  system  described  in  this 
report.  In  particular,  this  coordinate  system  has  one  axis  constrained  to 
lie  in  a  given  plane  and,  consequently,  does  not  rotate  with  the  missile. 

To  provide  a  working  simulation  and  to  determine  whether  this  new  co¬ 
ordinate  system  alleviates  the  conditions  producing  unsatisfactory  trajec¬ 
tories,  the  equations  derived  in  this  report  were  programmed  by  the 
Digital  Applications  Unit  for  the  IBM  709  computer,  Reference  5  contains 
a  description  of  the  corresponding  computer  program. 
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The  equations  again  consider  both  flat  and  spherical,  rotating  and  non-  i 

rotating  earth  cases,  and  again  make  use  of  the  Euler  transformations  *  ' 

(as  opposed  to  direction  cosines)  to  express  vector  quantities  in  various 
coordinate  systems. 

As  before,  guidance  factors,  motion  of  the  earth  along  its  orbit,  launcher 
effects,  and  asymmetric  missiles  are  not  included  in  the  derivations. 

These  represent  potential  areas  of  extension  of  the  present  equations  of 
motion. 


NOTATION 

Unit  vectors  in  each  of  the  three  orthogonal  directions  are  represented 
by  ij  ,  jjt  k;  where  subscript  i  denote^  the  coordinate  system  under  con¬ 
sideration.  Components  of  vectors  in  each  of  these  directions  will  have 
two  subscripts.  The  first  subscript  (X,  Y,  or  Z)  denotes  a  component 
along  the  i,  j,  k  axis,  respectively;  the  second  subscript  (I,  E,  etc.) 
denotes  the  referencing  coordinate  system.  Thus  VXE  is  the  component 

of  the  vector  V  along  the  i  axis  of  the  E-coordinate  system. 

Arrows  over  vector  quantities  denote  vectors  of  arbitrary  magnitude, 
bars  over  vector  quantities  denote  vectors  of  unit  magnitude. 

Subscripted  vectors  other  than  the  subscripts  mentioned  above  are 
enclosed  in  parentheses,  i.e.  (<3^). 


Finally  _ 


denotes  time  differentiation  of  vectors  relative  to  the 


.  th 
i 


dt 

coordinate  system,  while  a  dot  over  a  given  variable  denotes  scalar  dif¬ 
ferentiation  with  respect  to  time. 


3 


SYMBOLS 


Vector  Quantities 


‘i*  >i*  ki 
‘E’  >E*  kE 

iw>  kW 

:tf  !H’  kH 

*v>  i  V  k  v 
*M’  'm1  kM 

u 

» 

a 


oj 

Q 

(wT) 

-e 

CJ* 

-* 

R 

e 

V 

(Vw) 

If 

li 

•e 

J 

9 


(V) 


coordinate  systems 


Angular  velocity  of  (H)  relative  to  (I) 

Angular  velocity  of  (E)  relative  to  (I),  i.e.,  angular  velocity  of 
the  earth  about  its  axis 

Angular  velocity  of  (H)  relative  to  (E) 

Angular  velocity  of  (M)  relative  to  (H) 

Angular  velocity  of  (M)  relative  to  (I) 

Angular  velocity  of  (V)  relative  to  (E) 

Vector  from  center  of  earth  to  the  current  missile  CG  position 

-a 

Time  rate  of  change  of  R.  i.e.,  missile  velocity 
Vector  describing  wind 
Summation  of  forces  acting  on  missile 
Summation  of  moments  acting  on  missile 
Angular  momentum  of  missile 

Vector  giving  direction  of  gravitational  force  exerted  on  missile 
Velocity  of  missile  relative  to  the  air 
Thrust  mil  alignment  vector 


Thrust  vector 
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SYMBOLS  (Cent) 


Scalar  Quantities 


o,  o,  e,  v 

<jxj)' 

p 

d 

kr>A’  kN’  kF 
kH'  k<^p’  kd> 

rc 

an 

8 

Re»rth 

*o 

h 

*r  *a 
k 

a,  /I 

*vwl 

d> 

A,  B,  G,  H 


mass  of  rocket 

longitudinal  moment  of  inertia 
transverse  moment  of  inertia 
Eulerian  angles 

Angles  relating  (H)  to  (V)  coordinates 
Density  of  air 
Diameter  of  missile 
Aerodynamic  coefficients 

Aerodynamic  coefficients 

Distance  from  nose  of  missile  to  CG 

Distance  from  nose  of  missile  to  normal  center  of  pressure 
Distance  from  nose  of  missile  to  magnus  center  of  pressure 

Gravitational  acceleration 
Radius  of  earth 

Gravitational  acceleration  at  sea  level 

Altitude  of  rocket 

Thrust  malalignment  angles 

Portion  of  thrust  devoted  to  produce  jet  tortjue 

Anglos  defining  missile  position 

Magnitude  of  wind  velocity 

Angle  defining  direction  of  wind 

Angles  defining  initial  conditions 


SYMBOLS  (Cont) 


.  Denotes  time  differentiation 

"  Denotes  quantities  referred  to  (V)  coordinates 

>;  Denotes  the  vector  cross  product 


PROCEDURE 


Several  coordinate  systems  are  utilized  in  deriving  the  equations  of 
motion.  They  are  tabulated  here  for  later  convenience. 

i |,  jj,  kj  Inertial  coordinate  system;  origin  (O)  at  the  center  of  the 
earth,  k(  axis  coincident  with  positive  spin  axis  of  the 

earth,  ij  axis  coincident  with  0°  longitude,  and  jj  axis  so 
directed  as  to  form  a  right-handed  coordinate  system. 

i ,  j e ,  kE  Earth-Hxed  coordinate  sys  em;  identical  position  as 

ij,  jj,  kj  at  time  equal  zero,  but  coordinates  are  to  be  fixed 
to  the  earth. 


*W’  iw>  Azimuthal  coordinate  system;  this  coordinate  system  is  used 
to  introduce  wind  data  into  the  equations  of  motion.  The 
origin  of  these  coordinates  is  again  at  O,  the  iv  axis  is 

directed  towards  the  missile  CG,  the  k\y  axis  is  directed 
towards  the  positive  spin  axis  of  the  earth,  and  jw  is  so 
directed  as  to  produce  a  right-handed  coordinate  system. 

*H’  Ih*  Fixed  horizontal  plane  coordinate  system;  iH  axis  along 

missile  longitudinal  axis  is  directed  from  the  centroid 
towards  the  nose  of  the  missile,  is  constrained  to  lie  in 

a  horizontal  plane,  parallel  to  the  7F,  plane,  and  k~H  is 

directed  as  to  produce  a  right-handed  coordinate  system. 
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7  l 


.  n 

Fixed  vertical  plane  coordinate  system;  the  iy  axis  is  coin-  \| 

cident  with  the  i^  axis,  the  jy  axis  is  constrained  to  lie  per-  v 

pendicular  to  the  i£  axis  (i.e.,  remain  in  a  specified  vertical  * 

plane),  and  the  kv  axis  is  so  directed  as  to  determine  a  right- 
handed  coordinate  system. 

Missile-fixed  coordinate  system;  the  iy  axis  is  coincident 
with  missile  longitudinal  axis,  the  and  kM  axes  are  rigidly 
attached  to  the  missile  to  form  a  right-handed  coordinate 
system. 

The  procedure  will  be  divided  into  several  sections,  as  outlined  in  the 
Table  of  Contents,  each  describing  a  particular  aspect  of  the  equations  of 
motion,  with  the  last  section  combining  all  that  precedes  it  and  serving 
as  a  summary  of  the  equations  of  motion.  The  first  section  will  derive  the 
general  equations  of  motion,  and  will  include  the  necessary  transforma¬ 
tions  between  coordinate  systems. 

The  General  Equations  of  Motion 

As  usual,  Equations  1  and  2,  the  basis  for  Newtonian  mechanics,  pro¬ 
vide  the  foundation  for  the  equations  of  motion.1 


*V  >  i  V »  k  v 


‘M.  Im.  k  M 


-♦ 

If 


dtJ 


(i) 


dt 


(2) 


Strictly  speaking  Kquation  1  should  read: 


a 

VF 


A.  I  mV) 
dt 


d.V 

m  _! _  *  urn 


dt 


d’.R  -  . 

m  _ ! _  i  um 

dt*' 


<!') 


The  latter  term,  un,  involves  the  amount  of  mass  being  expelled  from  the  rocket  system, 
and  the  relative  velocities  of  the  rocket  and  the  exist  gases,  (u>.  This  term  is  known  as 
the  jet  reaction  and  is  considered  aspart  of  the  thrust,  whose  total  makeup  also  includes 

considerations  of  nozzle  deaign,  operating  temperatures,  and  pressures.  Kquation  l‘  uses 
u  in  place  of  the  seemingly  indicated  V  becajse  the  rocket  system  itself  hasnot  really 
been  defined  in  this  presentation.  The  interested  reader  is  referred  to  References  2  and 
}  for  a  more  complete  discussion. 
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Here  2F  and  2L  are  summations  of  forces  and  moments  acting  on  the  mis¬ 
sile,  R  is  a  vector  from  the  center  of  the  earth  to  the  current  CG  (center  of 

gravity)  position  of  the  missile,  and  J  is  the  total  angular  momentum  of  the 
rocket. 

Since  Equations  1  and  2  are  vector  equations  in  three-dimensional 
space,  three  component  or  scalar  equations  are  implicit  to  each.  To  deter¬ 
mine  explicit  directions  relative  to  the  missile  along  which  forces  and 
moments  can  be  conveniently  summed,  consider  a  coordinate  system  whose 
origin  is  at  the  missile  CG  and  whose  axes  are  directed  as  given  previously 

by  i,  j  k  _  Let  d  be  the  angular  velocity  of  this  coordinate  system 
H  H  n 

relative  to  inertial  coordinates,  where 


^  '  “xh‘h  +  “YH’H  +  UZHkH  ^ 

in  terms  of  the  (H)  coordinates. 

Utilizing  the  derivation  presented  in  Appendix  A,  Equation  1  may  be 
rewritten  as  follows: 


where  all  quantities  in  the  right  hand  side  of  the  equation  are  understood 
to  be  expressed  in  (H)  coordinates.  One  may  correctly  surmise  that  the 
presence  of  the  «  v  term  relates  the  motion  of  the  coordinate  system  (H), 
to  which  forces  and  moments  are  referred,  to  an  inertial  coordinate  system 
(I)  as  required  by  Newton’s  Laws  of  Motion. 

Rather  than  differentiate  Equation  4  directly,  it  is  convenient  to  let 


*  VXHiH  <  VYHjH  t  V7HkH 
dt 

so  that  Equation  4  becomes 


It 


dt 


m 


(V 

\  dt 


(5) 


(6) 
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Since  V  is  already  expressed  in  (H)  coordinates,  (see  Equation  5),  one 
can  write 

dHV  .  -  .  -  .  — 

— —  “  VXH*H  +  VYH>H  f  VZHkH  W 

dt 

Combining  Equations  3,  5,  6,  and  7,  and  performing  the  indicated  opera¬ 
tions,  one  then  obtains  the  three  component  force  equations,  namely: 

=  m[VXH  +  (“YHVZH  -«ZHVYH)1  in  the  rH  direction 
=  m[VYH  +(WZHVXH  "  WXHVZH)]  in  th<?  difeCtion  (8) 

=  +  ^WX H^YH  ~  “  YH^ XH^  in  the  kH  direction 

One  can  treat  the  moment  equation  (Equation  2)  in  a  similar  manner, 
obtaining  first 


IF 


XH 


IF 


YH 


IF 


ZH 


XL  =  ^  =  iil  r  cS  X  f  (9) 

dt  dt 

Denoting  (wx)  as  the  total  angular  velocity  of  the  missile  relative  to 

inertial  coordinates,  we  can  write  the  angular  momentum  in  the  following 
form: 


J  ^X(<JT}XH'h  +  IY(wT)YHiH  f  IY(toT)ZHkH 


(10) 


where  Ix,  Iy  are  the  moments  of  inertia  of  the  missile  about  the  longi¬ 
tudinal  and  transverse  axes,  respectively.  Note  that  the  assumption  of 
rotational  symmetry  (Iz  =  I  y)  is  implicitly  made  in  Equation  10. 

Since  J  is  already  expressed  in  the  (H)  coordinate  system,  the  time  rate 
of  J  can  be  expressed  as  follows:' 


^X^  T^XH*  H  '  ly(wx)YniH 

dt 


Iy(o.T)zH,tH 


(ID 


Again,  strictly  speaking,  terms  of  the  folloatng  form  should  be  included  in  the  differ* 
enti.ition  during  the  burning  stages: 

'^^Vxir  *y(wt)yh’  'y^'t^zh 

Howevei,  for  a  fairly  stable  rocket  (no  tumbling),  ly(  T^YIf  10 ^  ^Y*°',1^  ZH  c,n 
neglected.  The  expression  Iy(o)X)  Xj|  can  be  written  as  follows: 

ix(w-^XH  ink^or-^xH  *  m(0>T>XH  <k,)  2 

where  k  is  the  radius  of  gyration.  The  latter  term  can  be  neglected  because  usually  the 
burning  fuel  has  Tittle  effect  on  k\  In  addition,  for  low-spinning  rockets,  the  former  term 
also  can  be  neglected.  For  high-spin  rockets,  however,  it  k  <">x>x„  «*  combined  with 

other  terms  which  comprise  the  jet  torque  which  acts  to  produce  high- spin  rates.  The 
reader  is  referred  to  Reference  }  for  further  details  on  this  matter. 
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Appropriately  combining  Equations  3,  9,  10,  and  11  now  results  in  the 
component  equations  for  the  moments: 

^LXH  =  ‘^T^XH  +  (“yhM^T^ZH  ~  "ZhM^T^Yh) 

2L™  =  Iy(<y  t)  yh  +  ("zh^t)  xh  -  "xhM^tW  (12) 

SLZH  =  M^ZH  +  <"xhMw^YH  "  "YH^tW 

Equations  8  and  12  comprise  the  basic  six-degree  equations  of  motion. 

These  equations  are  of  limited  use,  however,  until  one  knows  explicitly 
how  the  (H)  coordinate  system  is  moving  with  respect  to  inertial  or  earth- 
fixed  coordinates.  Further,  one  often  has  vector  quantities  expressed  in 
earth  coordinates  (such  as  wind  and  gravitational  attraction)  which  must 
be  properly  introduced  into  the  (H)  coordinates.  Consequently,  additional 
relationships  between  the  (H)  and  (E)  coordinates  must  be  derived.  In 
particular,  one  must  know  how  the  (H)  coordinate  system  is  oriented  rela¬ 
tive  to  the  (E)  coordinates,  at  all  times.  The  orientation  of  (E)  relative  to 
(I)  must  also  be  known.  Although  other  methods  exist,  as  the  Introduc¬ 
tion  indicates,  the  use  of  Euler  angles  appears  to  be  the  most  straightfor¬ 
ward  for  the  present  study.  In  this  approach,  one  rotates  a  coordinate  sys¬ 
tem,  initially  coincident  with  the  (E)  system,  about  selected  axes  so  that 
after  the  rotations  are  performed  this  coordinate  system  will  have  the  same 
orientation  as  the  (H)  system.  One  will  then  have  the  means  for  expressing 
vectors  of  one  coordinate  system  in  terms  of  the  other. 

The  reader  is  reminded  that,  in  performing  these  rotations,  the  iH  axis 
is  to  be  coincident  with  the  missile  axis,  while  j ^  remains  parallel  to  the 

horizontal  i£  -  jp  plane. 

To  accomplish  this  end,  assume  the  existence  of  an  arbitrary  vector, 
extending  from  the  origin  of  the  (E)  system,  to  represent  the  missile  axis. 
First  rotate  the  (E)  system  about  the  kE  axis  by  an  angle  4>  so  that  the 
i'  axis  (ip  rotated)  coincides  with  the  projection  of  the  missile  axis  on 

the  i E  — j K  plane,  as  shown  in  Figure  1  (p  11). 
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Missile  axis 


Horizontal  plane 


Figure  I 

One  may  conveniently  express  this  rotation  in  matrix  form,  as  shown  in 
Appendix  B. 

i'  cos  0  sin  0  0  i 

j'  -  -sin  0  cos  0  0  )  (13) 


Another  rotation  is  yet  required  about  the  positive  j'  axis,  (j  E  rotated)  so 
tha:  iH  (i'  rotated)  is  coincident  with  the  missile  axis.  Denote  the  magni¬ 
tude  of  this  rotation  by  0,  leading  to  the  matrix: 

ijj  cos  0  0  -sin  0  i' 

rH  ■=  o  i  o  r  (i4) 

kH  sin  0  0  cos  0  k’ 

To  avoid  difficulties  later  on  resulting  from  ambiguity  in  the  determination 
of  initial  conditions  and  in  other  coordinate  system  transformations,  let  0 
be  defined  from  0°  <_0  <  360°  and  -90°  <_0  <.90°,  for  all  time!  Obviously 

special  attention  must  be  given  to  the  equations  of  motion  when  0  passes 
through  in/2  since  sharp  discontinuities  in  0  will  result.  Each  rotation 
matrix  is  a  linear  transformation;  hence,  one  may  obtain  the  (H)  coordinates 
directly  in  terms  of  the  (E)  coordinates  by  combining  both  rotations.  This 
is  equivalent  to  the  following  matrix  multiplication. 


‘As  an  example.  Figure  1  would  require  6  to  take  on  a  negative  value. 
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(15) 


- 

- 

“ 

~ 

~ 

~ 

'H 

cos 

e 

0 

-sin  6 

cos  p 

si  n  p 

0 

‘E 

4 

0 

1 

0 

-sin  p 

cos  p 

0 

i  e 

sin 

6 

0 

cos  6 

0 

0 

1 

kE 

producing 


*H 

|h 

kH 


cos  6  cos  p 
-sin  tp 
sin  Q  cos  ip 


cos  6  sin  \p  -sin  6 
cos  p  0 

sin  6  sin  p  cos  6 


*E 

>E 


(16) 


Note  from  Figure  1  (p  11)  that  jH  ( j '  rotated)  remains  in  a  horizontal 
plane,  or  is  perpendicular  to  kE.  This  is  reaffirmed  in  Equation  16  where 
the  component  of  jH  on  kE  is  seen  to  be  zero. 


One  may  obtain  the  inverse  of  Equation  16  to  obtain  (E)  in  terms  of 
(H)  as  follows: 


■ 

- 

• 

1 E 

cecp  -sp  seep 

4 

i  e 

ce  sp  cp  se  sp 

47 

kl. 

-so  o  ce 

-  kH 

(17) 


For  brevity  cos  A  and  sin  A  have  been  replaced  by  CA  and  SA,  respec¬ 
tively.  It  should  be  noted  that  Equations  16  and  17  actually  represent  three 
scalar  equations.  To  illustrate  this,  let  us  assume  the  velocity  components 
of  the  missile  expressed  in  the  (H)  system  are  known  and  we  wish  to  ex¬ 
press  the  velocity  in  the  (E)  system.  We  have 

^  "  VXH*H  *  VYH>H  4  VZHkH  "  VXE‘k  *  VYE*E  *  VZEkE  ( 18^ 


First,  obtain  i„,  j„,  kH  in  terms  of  iE,  jE,  kE,  which  can  be  done  by  using 
Equation  16.  Thus, 

i>(  -  C0C0iE  +  CB S0ji;  -  S0kE 

rn  -S0.7;  .  c0rE  *  ok;.  d9) 

kj,  soc^Tj  ♦  *  aik~ 
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It  then  remains  to  substitute  these  equations  into  Equation  18  and  equate 
the  coefficients  of  iE,  jE,  kE  to  VXE,  VYe.  VZE>  respectively.  Performing 
these  operations  will  yield: 


VXE  =  VxH0^'/'  "  VYHS^  +  VZH S0C0 
VYE  =  VXHC^  +  VYlf^  +  VzH56^ 


V 


ZE 


-vx„s«  .  vZHc« 


(20) 


This  is  the  required  transformation. 

Although  much  has  been  accomplished  with  Equations  13  through  20,  one 
must  go  further  to  obtain  the  rates  of  change  of  the  angles  0  and  6  as  the 
missile  position  varies  in  time. 

Let  Z  be  the  angular  velocity  of  (H)  relative  to  (E)  where,  in  (H)  coor¬ 
dinates, 

w  =  “XlfH  +  ^YH^H  +  k*  Ztfc  H  ^ 


Appendix  B  shows  Z  can  also  be  written  as 

Z  -  0k  E  +  Q\'  (22) 

Using  the  matrices  already  developed,  one  proceeds  by  expressing  both 
kE  and  j*  of  Equation  22  in  terms  of  the  (H)  coordinates  and  equating  the 
resulting  coefficients  of  iH ,  jH,  krt  to  the  coefficients  of  iH,  jH,  kH  in 
Equation  21.  This  will  produce 

~  -0  Sd 

WyH  =  b  (23) 

o»zH  ’  0C0 


The  rates  of  change  of  the  Euler  angles  become  simply 


WXH  a,ZH 

~W  ~ce 


(24) 
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One  may  note  that  several  angular  velocities  have  been  introduced  into 
the  equations  of  motion,  namely,  <3,  w,  and  (c3x),  1°  component  form,  this 

amounts  to  nine  unknowns,  a  completely  overwhelming  assignment  for  the 
three  moment  equations.  Fortunately,  not  all  these  quantities  are  inde¬ 
pendent.  It  can  be  shown  (see  Ref  4,  for  example)  that  if  a  coordinate  sys¬ 
tem  (H)  is  rotating  with  angular  velocity  dl  relative  to  a  coordinate  system 
(E),  and  (E)  is  rotating  with  angular  velocity  Q  relative  to  a  coordinate 
system  (I),  then  the  angular  velocity  of  (H)  relative  to  (I),  is  simply  the 
vector  sum  of  the  individual  angular  velocities.  Thus, 

w  =  Z  +  Q  (26) 


This  can  be  generalized  to  more  than  three  coordinate  systems.  For 
example,  continuing  to  use  the  definitions  given  in  the  Table  of  Symbols, 
we  are  also  at  liberty  to  write 

(eu-j)  =  Q  +  6=  n  +  5  +  Q  (27) 

It  is  now  necessary  to  obtain  explicit  expressions  for  3,  (wT),  and  (d>T), 
each  of  which  is  necessary  to  the  solution  of  the  general  equations  of 
motion. 

To  obtain  components  of  Q  in  the  (H)  system,  one  need  only  convert  kE 
into  the  (H)  coordinates,  since  the  kE  axis  was  assigned  to  be  coincident 
with  the  spin  axis  of  the  earth.  Therefore 

Q  =  Qk~  =  -QS0i~  +  QC0k^  (28) 


The  total  angular  velocity,  (wT),  of  the  missile  is  introduced  by  first 
making  use  of  the  (M)  coordinate  system.  Since  the  7^  and  if|  axes  are  to 
remain  coincident  at  all  times,  the  angular  velocity  fl  of  (M)  relative  to 
(H)  is  simply; 


n 


n  r 

XM  M 


n  r 

XH  H 


(29) 


By  the  use  of  Equation  27,  the  total  angular  velocity  can  now  be  written 
in  component  form  as  follows: 


(oiT)XH 

nXH  * 

WX  H  -  QS0 

(wT)  Yh 

WYH 

(30) 

Cti  _  . .  4 

QC0 

Zll 

14 


One  may  still  sense  the  existence  of  too  many  unknowns,  namely,  0,  and 
the  three  components  of  to.  Fortunately  Equation  23  gives  a  relation  between 
two  components  of  <o.  Eliminating  < ft  produces 

gj x h  ~  ~a>ZH  ^  ^  (31) 

Finally,  to  obtain  explicit  expressions  for  (d>T)  as  given  in  the  basic  mo¬ 
ment  equation,  Equation  31  is  substituted  into  Equation  30  and  the  latter 
is  differentiated  with  respect  to  time.  This  leads  to: 

(<VXH  =  ^XH-6iZHtan0-6iZHseC^6JYH-QC0^YH 

=  0JYH  (32) 

^T^Z.H  =  0J7.H  “  2  S0WYH 

This  completes  the  discussion.  The  equations  needed  to  obtain  the 
actual  trajectory  are  summarized  beginning  on  page  51.  The  next  section 
considers  the  forces  and  moments  which  constitute  the  left”  side  of  the 
equations  of  motion  (Equations  8  and  12). 

Forces  and  Moments 

This  section  considers  the  forces  and  moments  acting  on  the  missile. 
Excluding  guidance,  these  forces  and  moments  carv  be  divided  into  three 
categories  as  follows: 


1.  Aerodynamic 

2.  Gravitational 

3.  Jet 

Since  the  forces  and  moments  are  essentially  as  presented  in  Reference 
1,  only  a  brief  account  of  each  will  be  given  here. 


Aerodynamic  Forces 


Three  aerodynamic  forces  are  considered.  A  brief  table  giving  magni 
tude  and  direction  of  each  is  given  on  page  16. 
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Fore* 
Axial  drag 

Normal  force 
Magnus  force 


Magnitude' 

P*<V'kDA 
PW)2  kN 

pd,(a,T)Xf/Vf)k 


Di  faction 

Along  missile  axis 

Perpendicular  to  missile  axis 
in  plane  of  yaw 

Perpendicular  to  plane  of  yawJ 


*(V  )  is  the  magnitude  of  the  missile  velocity  relative  to  the  air. 

2The  yaw  plane  is  the  plane  that  contains  both  the  missile  axis  and 
resultant  velocity  vector. 


These  forces  are  illustrated  in  Figure  2. 


Figure  2 

The  axial  drag  acts  along  the  negative  i||  axis,  which  by  definition  is 
directed  from  the  CG  towards  the  nose  of  the  missile.  Explicitly. 

(Axial  drag)Xj|  /hF  ( Yf )'  k|M  (M) 

The  normal  force  components  act  opposite  to  the  directions  of  (V  )yM 
and  (V  ).  .  Schematically  we  have 


1(> 


Knowing  the  scalar  magnitude  of  the  resultant  normal  force  and  using 
the  geometry  of  Figure  3,  one  can  deduce  the  component  of  this  force  act¬ 
ing  along  jH,  namely 


(Normal  force)  YH 


-pdJ(Vr)JkN  cos  a 


-pd’(Vr)1  kN(Vr)yH 

V(v  ry  +  (vry 

v  'yh  Z  H 


(34) 


Similarly,  in  the  kt!  direction 


-pd3(Vr)JkN(Vr)ZH 


(Normal  Force)7,, 


(35) 


The  Magnus  force  acts  in  a  direction  perpendicular  to  the  plane  of 
yaw,  as  indicated  in  Figure  4.  It  should  be  noted  that,  for  reversed  spin, 
the  Magnus  force  will  act  in  the  opposite  direction. 


Using  the  scalar  magnitude  and  the  geometry  of  Figure  4  above,  it  is 
readily  ascertained  that  ir.  the  jH  direction 


~pd5 ((/.>_) Y..(V  )kP(V  )7u 

(Magnus  force)Y}|  -  ~pd3(<oT)XH(V  )kf  siri  a  - I_±ll — i — _ — i — —  (36) 


nWyh  3  <v»>zn 


Similarly,  in  the  k„  direction 


(Magnus  force) 


pd1((uT)X||(V  )kF(V)Yn 


V(v  y  (  (v  y 

'  rill  (711 


(37) 


Aerodynamic  Munivnf i 

Proceeding  in  a  similar  vein  as  for  the  aerodynamic  forces,  the  following 
table  gives  the  magnitude  and  direction  of  each  of  the  aerodynamic  mo¬ 
ments  considered: 


ZH 
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Moment 

Magnitude 

Direction 

Restoring 

pdJ(V)JkN(rc  -  An) 

Perpendicular  to  the 
plane  of  yaw 

M  agnu  s 

Pd>("T)xH(Vr)kF(rC  "  AM> 

In  plane  of  yaw 

Yaw  damping 

pd4(V)  \/(w/YH  +  (<JT)i,KkH 

Perpendi  cular  to 
missile  axis 

Roll  damping1 

pd4(Vf  )k^p(^T)  XH 

Along  missile  axis 

Roll2 

Along  missile  axis 

while  for  fin-stabilized  missiles  k^  is 

to  be  changed  to  k^. 

2 

The  roll  moment  is  used  only  for  missiles  with  canted  fins. 


‘V 


s  used  for  spin-stabilized  missiles, 


In  obtaining  the  components  of  the  above  moments,  one  may  use  the 
relation  "distance  *  force  -  moment,’  ’  applied  in  a  manner  consistent 
with  the  geometry  of  the  system.  When  expressed  in  vector  form  the  restor¬ 
ing  moment  equation  becomes 

(Restoring  moment)  (rc  -  A H  *  [(normal  force)YHjH  +  (normal  force)ZHkHl  (38) 
This  produces,  in  component  form 


(rr  -AN)pd2(Vr)JlcN(Vr)ZH 
(Restoring  momentJyjj  — ^  »  '  — 

v'<V)‘v„  *  (V)'z„ 


(59) 


in  the  in  direction,  and 


(Restoring  moment)  ZJj 


(40) 


'<v,»‘yh  •  <v.>; 


7M 


in  the  k,,  direction.  Note  that  both  magnitude  and  direction  are  consistent 
with  the  above  table. 

Doing  likewise  for  the  Magnus  moment  produces  along  the  j|{  axis 


(Magnus  moment)Ytl 


“(rC  ~  AM)^»d 1  ( oj  T)  x vr  vr  )y  H 


(41) 
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while  along  the  kH  axis 


(Magnus  moment) 


<*r 


AM)pd3(a|T)xH(V  )kF(V  )ZH 
Vn)yH  ;  ^>ZH 


(42) 


The  yaw  damping  moment  acts  to  diminish  the  yaw  of  the  missile; 
hence,  this  moment  acts  opposite  to  the  lateral  angular  velocity  of  the 
missile.  The  components  of  this  moment  are  readily  ascertained  to  be 

(Y aw  damping)  Y^  -  -pd4(V)  cot)yh  ^  ^T^ZH  ^  Y  (43) 

(Yaw  damping)  ZH  =  ~pd4(Vr )  \/(wT)YH  +  (wT)zH  kH  sin  8y  (44) 

where  5Y  is  shown  in  Figure  5  below. 


Resultant  lateral 
angular  velocity 

\Z(«UT)YH  +  (wt^ZH 


Writing  cos  SY  in  terms  of  U>T)YH  and  (wY)ZH 


Cos  SY  _  (<"T>YH 

v("JY)Yh  1  (wj)zu 


(45) 


and  substituting  into  Equation  43  leads  to  the  component  along  i,. 


(Yaw  damping)yH  ~pd*(V )  (wT)YHk(1  (46) 

similarly,  for  the  component  along  FJ(, 

(Yaw  damping), ()  -pd4 ( V  )  (wT)7Hlt(|  (47) 
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The  roll  damping  moment  acts  to  reduce  the  spin;  hence,  it  is  introduced 
into  the  equations  of  motion  with  a  minus  sign  prefixed. 


The  direction  of  the  roll  moment  depends  upon  the  orientation  of  the 
canted  fins.  If  the  cant  produces  positive  spin  (i.e.,  a  clockwise  rotation 
of  the  missile  looking  from  the  rear  of  the  rocket),  then  this  term  is  intro¬ 
duced  with  a  plus  sign.  Like  the  roll  damping  moment,  the  direction  is 
along  the  longitudinal  axis  of  the  missile. 

It  should  be  noted  that,  since  none  of  the  aerodynamic  coefficients 
(kDA,  kN,  kF,  etc.)  are  assumed  to  be  linear  in  nature,  they  are  not  to  be 

taken  as  slopes  to  be  multiplied  by  angle  of  attack.  Rather,  these  coeffi¬ 
cients  are  point  values  obtained  directly  from  experiment  as  functions 
both  of  Mach  number  and  angle  of  attack. 

Table  1  (p  22)  summarizes  in  component  form  the  aerodynamic  forces 
and  moments  considered  thus  far. 

Gravitational  Fore* 

In  evaluating  this  force,  distinction  must  be  made  between  spherical 
and  flat  earth  cases.  In  the  latter  case,  the  gravitational  force  acts  in  a 
constant  direction  (independent  of  the  missile  position)  while  in  the  for¬ 
mer  case,  the  gravitational  attraction  is  directed  towards  the  earth’s  cen¬ 
ter  from  the  missile  CG.  In  both  cases,  specific  expressions  can  be  de¬ 
rived;  the  simpler  case  will  be  treated  first. 

Flat  Earth  Ca*o 

In  the  absence  of  the  earth’s  rotation,1  both  theT,  and  kE  axes  are 

taken  as  extending  in  positive  vertical  direction.  The  origins  of  the  (I) 
and  the  (E)  coordinate  systems  are  both  located  on  the  surface  of  the 
’’flat”  earth.  Hence,  mgy  (y  being  a  unit  vector  representing  the  direction 
of  the  gravitational  force  acting  on  the  missile)  becomes. 

mgy  -mgkj  -mgkj.- 


tine  should  have  little  desire  to  reline  the  model  to  include  the  earth  s  rotation  and 
yet  allow  the  assumption  of  a  plane  to  represent  the  earth  s  jjeometr>. 
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Replacing  kg  by  its  representation  in  (H)  coordinates,  similar  to  Equa¬ 
tion  28,  one  obtains 

mgy  -  -  mgCOk^  (49) 

in  a  form  suitable  for  substitution  into  the  general  equations  of  motion. 

Spherical  Earth  Case 

For  this  case,  the  direction  of  the  gravitational  force  depends  upon  the 
missile  position.  In  particular,  let  the  current  missile  position  be  repre¬ 
sented  by  a  vector  R,  where 

R  =  RXe‘E  +  RYE>E  +  RZEkE  ^ 


as  shown  in  Figure  6. 


Dividing  R  by  its  magnitude  produces  the  unit  vector 

R  A-  RXE‘K  1  RYK'E  *  R7.Ek k  (51) 

I  ft  I  VRxg  *  Ryf  *  rze 

Clearly, 

mgy  -mgR  (52) 
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Combining  Equations  51  and  52  and  obtaining  i£,  j£,  k£  in  terms  of  the 
(H)  coordinates,  there  results: 


,  i 


(Gravity  force)xH  =  - - [R  +  R  ygC0S0  -  R  z 1 

vRXE  +  rye  +  RZE 


(Gravity  force)  yh  =  - 1212 - [ -Rxe^0  +  RyeC0 ]  (53) 

^RXE  +  RYE  +  RZE 

(Gravity  forceJ^H  =  - — — - [  R  XeS0C0  +  RyE^S0  +  R  ze£^1 

^RXE  +  RYE  +  RZE 


One  should  note  that  the  gravitational  acceleration,  g,  is  a  function  of  alti¬ 
tude.  In  both  cases,  one  may  write 


g 


Rearth 

h 


(54) 


where 


go  -  acceleration  of  gravity  at  sea  level, 
Rearth  =  ra<l‘us  °f  chc  earth, 

I  -  current  altitude  of  missile 
h  <  -  (R  ZE  ^or  ^at  earfh  case) 


+  RYE  +  RZE  for  the  spherical  earth  case. 


Since  the  gravitational  force  acts  at  the  CG  of  the  missile,  there  are  no 
moments  associated  with  this  force.  Table  2  (p  25)  summarizes  the  results 
of  this  derivation  in  component  form. 


J«t  Force* 


A  single  jet  force  is  considered,  the  thrust,  which  imparts  forward  motion 
to  the  missile.  Ideally,  this  force  should  act  along  the  missile’ s  longitudi¬ 
nal  axis;  however,  due  to  imperfections  in  the  rocke.  design,  the  actual 
nature  of  propellants  and  other  factors,  there  may  be  a  component  of  thrust 
perpendicular  to  the  longitudinal  axis.  Since  the  thrust  vector  is  defined  at 
a  given  time  relative  to  the  missile  body,  the  thrust  components  will  first 
be  .cified  in  the  (M)  coordinate  system  (i.e.,  the  coordinate  system  that 
is  rigidly  attached  to  the  missile). 
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Using  the  symbols  presented  in  Figure  7,  one  can  easily  ascertain  these 
components  to  be  as  in  Equation  55. 


Figure  7 


(Thrust)^  =  -  T  cos 

(Thrust)yM  =  ^YM  T  sin  5y  cos  <5A  (55) 

(Thrust)^  -  ^7M  -  T  sin  sin  5A 

It  is  necessary,  of  course,  to  determine  the  components  of  thrust  in  the 
(H)  coordinate  system  along  whose  axes  the  forces  are  summed. 

Figure  8  shows  a  typical  relationship  between  the  (H)  and  (M)  co¬ 
ordinates. 


Figure  8 
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Here  S  is  defined  as  the  angle  between  the  jM  and  j H  axes.  It  is  readily 
seen  that 


hi 


1  0 

0  cos  S 
0  -sin  S 


- 

- 

0 

sin  S 

cos  S 

kH  . 

and,  from  Equation  56, 


«M  =  *H 

j m  =  cos  S  )H  +  sin  SkH 
=  -sin  S  j^j  +  cos  S  kjj 


(56) 


(57) 


Combining  Equations  55  and  57  produces  the  required  result: 


T 

1  XH 

T  cos 

tyh  = 

T  sin  cos  ( 5^  +  S) 

(58) 

tzh  = 

T  sin  8  j  sin  (5  ^  +  S) 

Jet  Moment* 

Two  jet  moments  exist.  One  unintentionally  arises  from  the  fact  that  the 
line  of  thrust  may  not  pass  through  the  centroid  of  the  missile.  The  other 
is  an  intentional  jet  torque  which  causes  the  rocket  to  spin  about  its  axis 
of  revolution.  To  be  explicit  about  the  former,  one  may  define  a  vector? 
from  the  CG  of  the  rocket  to  the  point  of  application  of  the  thrust  at  the 
nozzle  exit.  As  before,  this  vector  is  defined  in  the  missile  coordinate  sys 
tern,  i.e., 

*  rXM'M  '  rYM»M  *  r7.MkM 

The  corresponding  jet  moment  is  now  given  by 

(Thrust  moment)  ■  r  v  (Thrust) 
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which  in  component  form  becomes: 


(Thrust  moment) 

rYMTZM  “  rZMTYM 

(Thrust  moment)  ym 

rZMTXM  ~  rXMTZM 

(Thrust  moment) 

rXMTYM  “  rYMTxM 

Converting  these  components  into  the  (H)  coordinates  and  substituting 
Equation  55  finally  produces 

(Thrust  moment)XH  -  T  sin  ^flryM  s*n  ^A  ~  rZM  cos  ^A^ 

(Thrust  momentiyn  T  cos  ^  T  ~  rXM  s*n  ^  T  s>n  8  cos  S 

-  T  [r  sin  cos  8^  -  rYM  cos  sin  S  (60) 

(Thrust  moment  >ZH  T  trZM  cos  -  rXM  sin  6.^.  sin  SA1  sin  S 

\  T  Ir^^  sin  8j  cos  8^  -  rYM  cos  Sjl  cos  S 

The  final  moment,  the  jet  torque,  will  be  assumed  proportional  to  the 
thrust.  If  this  moment  produces  positive  spin,  this  term  is  to  be  intro¬ 
duced  into  the  equations  of  motion  with  a  plus  sign  prefixed.  Mathemati¬ 
cally, 

(Jet  torque)  k  T  i  (61) 

This  concludes  the  discussion  of  forces  and  moments. 

One  might  note  that  V  ,  the  velocity  of  the  missile  relative  to  the  air, 
has  not  been  defined  in  the  presence  of  winds.  In  the  next  section  of  this 
report  explicit  expressions  will  be  derived  for  this  vector  in  the  (H)  coor¬ 
dinate  system. 

Introduction  of  Wind* 

This  section  derives  expressions  that  can  be  used  to  introduce  wind 
into  the  equations  of  motion  and  thus  determine  (Vf).  For  the  spherical 
earth  case,  it  will  be  assumed  that  wind  data  exists  at  preselected  points 
on  the  earth’s  surface.  For  a  flat  earth,  ii  is  assumed  that  wind  is  pre¬ 
sented  as  a  function  of  range  and  altitude. 

In  both  cases,  if  one  knows  the  velocity  of  the  missile  (VM)  and  the 
and  the  air  velocity  ( V A),  each  relative  to  the  earth's  surface,  then  the 
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velocity  of  the  missile  relative  to  the  air  (V  )  is  given  by 


<Va> 


(62) 


For  the  flat  non-rotating  earth,  wind  is  likely  to  be  given  relative  to  the  ground; 
therefore,  use  of  the  (E)  coordinate  system  is  appropriate  and  Equation  62  may 
be  written  as 


(V  )  (Vx,  -(V^)x|.)i{.  '  (VYK  -(V^)Yp)ii.-  *  (V/p  (Va)/.K)k  f.  (63) 
where  (V^)  is  the  actual  wind  vector. 


It  is  an  easy  matter  to  convert  Equation  63  into  the  (H)  coordinates. 
Doing  so  produces 


(V  )XH  I VX1,  -  (Vw)xi,  1  coco  *  lvYI;  -(VW)YK1(:OS0  -  |vZK  - <vw)ZEls/> 

( V  )  Y u  “(Vj-p  -  ( \'^’)xp  *  l^yp  (^w)  Y  v)f (64) 

(V,)/h  l  V  x  |.-  -  ( V^)xE  +  lVYp;  -  (V  \/j)  Yp]S0Si^  t  [VZp  -  (V^)Zp)C() 

Of  course,  (V«-)XK*  (V»-)yk.  and  (V«')zk  can  be  given  as  functions  of 
range. 

The  spherical  earth  case  is  a  bit  more  complicated  because  of  the 
geometry.  Also,  one  must  account  for  the  motion  of  the  air  induced  b_*  the 
earth’s  mtation  in  addition  to  wind  (i.e.,  disturbance  of  the  air  within 
this  rotating  air  mass). 

Recalling  that  Equation  *i 


d.,R 


.  w  •  H 


dt 


VXll'll  1  vyh'u  '  v / II *  II 


is  the  velocity  of  the  missile  relative  to  inertial  coordinates  (by  definition 
of  u),  one  must  also  express  the  motion  of  the  air  induced  b)  the  earth  s 
rotation  relative  to  inertial  coordinates.  By  so  doing  one  can  obtain  the 
.... i , .  c ... ; ; i..  ...  i. ...  »s«.  '■o'atinEair  mass .  Th  is  will  be  the 

equivalent  of  Equation  62  referred  to  inertial  coordinates. 


The  air  mass  rotates  as  a  whole  with  the  earth,  this  is  shown  physically 
by  noticing  that  the  rotation  of  the  air  mass  is  not  apparent  to  one  on  the 
earth's  surface.  The  velocity  of  this  rotating  air  mass  at  the  current  mis¬ 
sile  position  (specified  by  the  vector  R)  is  simplv 


The  velocity  of  the  missile  relative  to  the  rotating  air  mass  therefore 
becomes 


vM/Rot.  air  mass 


v  -  (3  *  ft) 


(66) 


One  must  now  account  for  the  wind  itself.  For  the  spherical  earth  case, 
winds  may  be  conveniently  introduced  into  the  equations  of  motion  by  use 
of  the  azimuthal  coordinate  system  (W).  This  coordinate  system  is  de¬ 
scribed  on  page  6  and  is  illustrated  in  Figure  9  (p30),  where  the  origin  of  (W) 
has  been  translated  to  the  earth’ s  surface  at  O'  for  visual  purposes.  The 

A 

wind  can  be  considered  to  be  a  function  of  the  height  h,  of  the  missile  CG 
above  the  earth’s  surface,  (see  Figure  9).  One  method  of  designating  the 
wind  in  the  (W)  coordinates  is  to  specify  wind  magnitude  and  wind  direction 
measured  from  North  at  given  altitudes.  It  will  be  assumed  in  this  report  that 
the  wind  vector,  V\(h),  has  no  vertical  component.  Thus  at  a  given  altitude 
Vw(fi)  shou.d  be  expressible  only  in  the  iw~ plane,  as  presented  in 
Equation  67,  using  the  geometry  of  Figure  9. 


V) 


y '  “W* 


(Vw)  -  |v^)|  sin  <£  j  y  +  |V^  (£)l  cos  4>  k\f/  =  (V^)y^j*jp  +(^w^ZW^W  (67) 

The  presence  of  ft  emphasizes  the  dependence  of  wind  on  altitude.  The 
point  O'  on  the  earth’s  surface  is  defined  by  a  and  y  where 


a 


y  --  tan 


f  R7.E 


(68) 


where  if 


rxe  0 


RXE  0 


ryk>0 

0 

<a  <90° 

rye  <0 

270° 

_  a  <  360 c 

R  Y  K  '■  0 

90° 

<  a  <  180 

ryk  0 

180° 

<  a  v  270 

and 


-90  ”  _  y  _90° 


It  remains  to  obtain  the  wind  components  i.i  the  (H)  coordinate  system 
and  add  the  results  to  Equation  66.  This  is  done  by  first  obtaining  jw 
and  in  terms  cf  (E)  and  then  (by  the  now  overworked  Equation  17) 
obtaining  (E)  in  terms  of  (H). 

Rotating  lirst  about  k{.  through  the  angle  «  and  then  rotating  about  j' 
through  the  angle  ft  {ft  -v  in  Fig  9).  one  can  see  that  the  resultant  co¬ 
ordinate  system  will  be  coincident  with  (W). 


Performing  the  indicated  operations  results  in 


CftCn  C.ftSn  -$ft 

*  K 

-Si <i  (■ a  0 

Sft(  it  Cft 

\  . 

31 


(69) 


Combining  Equation  69  with  Equation  17  completes  the  transformation. 
Thus , 


lw 

C/SCa  C/SSa  -Sj8 

" 

COCi/j  -S  0  S0C0 

j  w 

= 

-Sa  Ca  0 

C0S0  C0  S0S0 

kw 

S/3Ca  S/SSa  C/3 

-S(9  0  C<9 

kH 

leading  to 


r  1 

“ 

*w 

^11  ^12  A,j 

5 

i 

A.  A  .  A, , 

t 

Jw 

“2  1  “2  2  “2  J 

'H 

^Jl  ^12 

1 

u1 

where 

An  =  C/3CflC0C0  +  CpSaCdSi/,  +  S/3S0 

A13  -CfiCaSi/j  +  C/3S aCifi 

Als  •  C0CaS0Cdf  f  C.ftSaSOSi/j  -  S/3C0 

A,,  =■  -SaCOCi/j  t  CaC0S0 

Aj  j  SaSi/»  4  C aC0 

Aj5  ShS9Cu  .  CaSDSi.) 

A„  S/JCciC0C0  t  SfiSaCO  S0  -  C/3S.9 

A„  -SfiCaSiji  •  SfiSaC^j 

A„  SflCuSOCO  *  SflSaSflStf  .  C/K9 

Clearing  up  the  remaining  threads,  we  have 

(^w)  (V  w)  '  (Vv) ^  VW  ^XH*H  *  fVvi  Jyn'n  '  ^VW^/UlcH 
where  the  usual  algebraic  manipulation  results  in  the  necessary  equation 


(\7XII  »  Aji  ’  (VwV.w  A>> 
(VV..  (VV*  A»  ’  (Vzw  A» 
(V/.H  (Vy*A"’(V/..A« 


(72) 


(73) 


Finally,  the  sought  after  vector  is 

(Vr)  V-(QvR)-(Vw) 

* 

Note  again  that  Vv  prescribes  the  wind  velocity  relative  to  the  rotating 

air  mass.  Thus,  for  a  missile  moving  on  or  above  the  surface  of  the  earth 
with  a  velocity  of  Q  ^  $  (stationary  orbit  or  fixed  surface  point),  in  the  ab¬ 
sence  of  winds  would  cause  the  relative  velocity,  Vf  to  be  zero.  Writing 
Equation  73  in  component  form  yields 


(VxH 

^XH  (V\h')xh 

(Vyh 

vVH  aC^Rx,,  +  Q^R7h  -(Vw)yh 

(74) 

lVr>ZII 

V7.H  f  nS^RyH-  ^VzH 

;here,  of  course,  the  magnitude  of  (Vr)  is 


W  V<y>*xfi 

This  completes  the  analysis  of  winds. 

Initial  Condition* 

Initial  conditions  for  the  spherical  earth  case  will  be  derived  first. 
Initial  conditions  for  the  flat  earth  are  immediately  obtainable  upon  spe¬ 
cializing  certain  of  the  parameters. 

In  order  to  start  the  trajectory,  a  compete  set  of  initial  conditions  must 
be  provided.  These  are  tabulated  below  for  convenience. 


xt: 

VXH 

0(0) 

**XH 

v  l- 

VY». 

0(0) 

<|JYH 

V/H 

<■'/!! 

Those  conditions  that  depend  upon  a  coordinate  system  for  representa- 
non  can  bo  fiv«  in  rerms  of  «(.«  oordinarc  sysrems.  ,1  des.red,  as  Ion* 
as  known  transformation  equations  exist  between  such  coordinate  systems. 

Of  the  conditions  stated  above,  the  mtssil.  or.en.ation  expressed  in 
terms  of  )',<(»  and  »(#>  is  not  usual!)  known,  per  se.  But  the  follow, n* 
measurable  quantities  ate  (or  al  Ic"  r  'hould  bet  known: 


Longitude  of  the  launch  point 
Latitude  of  the  launch  point 


A° 


B° 

Azimuthal  heading  of  the  missile  G°  (This  will  be  defined  more  explicitly 

later  in  this  report.) 

Angle  of  declination  of  the  missile 
from  the  local  vertical  H° 

It  is  possible,  as  will  be  shown,  to  express  0(0)  and  0(0)  as  a  function 
of  the  variables  A°,  B°,  G°,  and  H°.  To  do  this,  one  rotates  a  coordinate 
system  initially  (t  0)  coincident  with  (E)  through  the  above  measured 
angles.  Upon  completion  of  the  rotations,  the  resultant  coordinate  system 
is  identical  in  orientation  with  the  (H)  coordinate  system  at  launch.  The 
(H)  system  is  also  defined  by  simply  specifying  0(0)  and  9(0).  Each  of 
the  resultant  rotations  can  be  expressed  in  matrix  form,  one  matrix  contain¬ 
ing  the  four  known  variables  and  the  other  matrix  containing  the  two  un¬ 
known  variables.  One  is  then  at  liberty  to  equate  corresponding  terms  of 
the  two  matrices,  since  each  actually  represents  the  same  coordinate  sys¬ 
tem.  A  particular  notational  scheme  is  used  to  distinguish  the  several 
rotations  performed.  Attached  to  each  unit  vector  is  a  series  of  primes,  the 
number  of  which  denotes  the  number  of  rotations  already  performed. 

The  purpose  of  the  first  two  rotations  is  to  locate  the  i  axis  so  that, 
when  extended,  it  passes  through  the  CG  of  the  missile  at  time  t  -  0. 

This  is  done  by  first  rotating  about  k through  the  angle  A°  so  that  the 

i'  axis  lies  on  the  same  longitude  as  the  missile  CG  and,  secondly,  by 
rotating  about  j'  through  the  angle  B",  the  latitude  of  the  missile  position. 
These  rotations  are  u^ed  in  the  same  order  to  obtain  the  position  of  the 
azimuthal  coordinates;  hence,  one  may  use  Equation  68  with  «  and  (3 
replaced  by  A  and  B,  respectively.  Thus 


The  purpoi-c  of  the  next  two  rotations  is  to  orient  the  i’1  axis  so  that  it 
becomes  coincident  with  the  missile  axis. 

In  particular,  the  third  rotation  is  about  i'1  through  the  angle  G",  so 
that  after  this  rotation  is  performed  -k’ ' '  coincides  with  the  projection  of 


-a 


the  missile  axis  on  the  j *  * '  —  k ’  1  *  plane.  In  matrix  form,  we  have 


The  fourth  rotation  is  about  j'  1  '  through  the  angle  H°  so  that  i' ' '  '  is 
finally  coincident  with  the  missile  axis.  Upon  examining  the  geometry  of 
these  rotations,  11°  is  seen  to  be  the  angle  of  declination  the  missile  makes 
with  the  local  vertical  i' ' . 

Mathematically, 


Combining  all  rotations  performed  results  in 


whe  re 


b„ 

CHCBCA 

-  SHSGSA  - 

SHCGSBCA 

b„ 

CM CBS A 

t  SHSGCA  - 

SHCGSBSA 

b„ 

-C  MSB  - 

SHCGCB 

b,. 

-cgsa  t 

SC.SBCA 

bj  a 

CXiCA  .  : 

SGSBSA 

bj  i 

SGCB 

b»i 

SIK  Me  A 

.  CIISGSA  ■ 

(  1K  X.SIU  A 

l*„ 

SM(  MSA 

(  llSt.t  A  • 

(  IK  C.SMSA 

b„ 

-SMSH  • 

(  IK  C.(  M 

Although  one  may  feel  that  a  sufficient  number  of  rotations  have  already 
been  performed  to  obtain  i *  *  *  '  coincident  with  the  missile  axis,  there  is 


yet  no  assurance  that  j""  will  lie  in  a  horizontal  plane,  as  required.  This 
situation  is  remedied  by  performing  a  fifth  rotation  about  i  ""  through  an 
angle  L  such  that  the  resultant  position  of  j""  lies  in  the  horizontal 
plane.  In  matrix  form 


r 

r- 

*H 

1  o  o 

'fill 

in 

= 

0  CL  SL 

j  MM 

0  -SL  CL 

k"" 

The  magnitude  of  L  is  as  yet  unknown.  Substituting  Equation  78  into 
Equation  79  produces  the  matrix 


b„ 

bj  ,CL  4  bj,SL 
-bj,SL  t  bs,CL 


bn 

b2  j  CL  4  b32  SL 
-b2  2  SL  f  b12  CL 


b2  jCL  4  bj,SL 
-b2  jSL  4  b,jCL 


(80) 


which  enables  one  to  convert  from  (H)  to(E)  coordinates.  Equation  16 


*H 

COCO  COSO  -SO 

1 E 

-so  co  o 

»E 

kn. 

S0C0  S0S0  co_ 

kF  . 

also  relates  the  ( II)  and  (E)  coordinates.  Since  these  sets  of  equations 
are  equivalent  at  the  start  of  the  trajectory,  one  simply  compares  corres 
ponding  terms,  from  which  one  obtains  L  as  well  as  0(0)  and  0(0). 

In  particular, 


Sm  0(0'  -b,, 


Sin  0(0)  - _li _ 

Cos  0(0) 


Cos  0(0)  - U - 

Cos  0(0) 


(81) 


(82) 


My 


and  by  using 


b23  cos  L  +  bn  sin  L  ■=  0 
bjj  cos  L  i  bjj  sin  L  =  C0 


(83) 


one  obtains  for  L 


Cos  L  =  b» 

bj  j  bJ2  -  bSJ  b3  j 

(84) 

Sin  L  =  b:a 

b2  a  bJ2  -  bn  b2  j 

Although  L  need  not  be  explicitly  used  in  determining  0(0),  0(0),  L  com¬ 
pletely  describes  the  missile  orientation  at  launch. 

Since  0  was  restricted  between  ±90°,  the  quadrants  of  each  angle  are 
uniquely  determined  in  the  preceding  equations. 

The  initial  linear  and  angular  velocities  may  be  given  in  terms  of  either 
the  (E)  or  the  (H)  coordinates,  although  w  is  usually  most  easily  expressed 
in  the  (H)  coordinate  system.  Perhaps  it  would  be  best  at  this  point  to 
illustrate  the  theory  by  a  brief  example. 

Let  us  suppose  that  it  is  desired  to  launch  a  rocket  as  pictured  in  Fig¬ 
ure  10  (p  38).  We  make  the  following  assignments: 

A  =  30°  (Longitude) 

B  -  -40°  (Latitude.  Note  the  negative  rotation  about  j'  for  the  Northern  Hemi¬ 
sphere;  hence,  the  minus  sign). 

G  =  10°  (190°  —  180°  so  that  the  negative  k'"  axis  will  coincide  with  the 
line  0*  P). 

H  --  35°  (Angle  of  declination  from  the  local  vertical). 

0(0)  and  0(0)  can  then  be  determined  by  use  of  the  equations  up  to  and 
including  84. 

Let  us  now  introduce  a  5°  yaw  angle  in  a  plane  containing  the  local 
vertical,  measured  towards  the  north,  as  pictured  in  Figure  11  (p  39).  From 
the  geometry  of  Figure  11  (b),  one  can  obtain  the  components  of  the  mis¬ 
sile  velocity  in  the  (E)  system  and  use  the  transformation  Equation  17  to 
obtain  V  ,  V  ,  and  V  .  fi  is  simply  the  spin  of  the  missile  as 

All  iM 
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Local  vertical  ( i  "  ) 


Diagram  Illustrating  Initial  Conditions  (II) 


North 


Local  vertical 


Initial  velocity  vector 
Missile  axis 


Local  vertical 


Velocity  vector 
Missile  axis 


Figure  1 1 

Y) 


observed  in  the  CH)  coordinate  system.  and  wZH  are  the  lateral 
angular  velocities  of  the  missile  as  observed  from  the  surface  of  the  earth. 
Sketching  the  rotations  in  sequence  before  numerical  values  are  assigned 
to  these  quantities  can  be  helpful. 

We  will  now  specialize  the  results  so  that  they  are  appropriate  to  the 
flat  earth  case.  As  has  been  mentioned  previously,  the  vertical  is  to  be 
the  kE  axis,  while  the  ip  axis  will  be  taken  as  downrange.  One  conveni¬ 
ent  way  of  ascertaining  initial  conditions  is  to  first  obtain  the  i"  axis 
coincident  with  the  kE  axis;  then,  with  Gr  and  H°  known,  compute  yi(0) 
and  0(0).  Referring  to  Figure  10  (p  38).  we  replace  the  30°  with  0°  and 
the  -40°  with  -90°  (or  270°).  The  i"  axis  thus  achieves  its  first  desired 
orientation.  Mathematically, 

A  0° 

B  270° 

Our  matrix  Equation  78  for  arbitrary  angles  G  and  H  now  becomes 


["  i"" 

SHCG 

SHSG 

CH 

‘E 

-sc. 

CG 

0 

!e 

(85) 

k~ 

-CHCG 

•OiSG 

SH 

.  kE  - 

When  the  appropriate  quantities  are  substituted  imo  Equation  83,  L  be¬ 
comes  zero.  Using  these  results  to  obtain  <£{0)  and  0(0)  for  the  flat  earth 
case  gives  the  following: 


u(0)  G 

0(0)  H  -  90  (86) 

L  O'* 

for 

A  O' 

b  :?o 

The  orher  quantities  VX(l . '‘V.H  arc»  as  before,  treated  analogously, 

although  R  x  R  Y  p,  and  R  / ,  may  (if  appropriate)  be  assigned  values  of 

zero. 
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This  concludes  the  section  on  initial  conditions.  More  information  will 
be  provided  on  '  is  subject  when  the  alternate  set  of  equations,  as  derived 
in  the  next  sec.  on,  is  used. 

Singularity  Conditions  ana  an  Altornato  Sot  of  Equations 

This  section  presents  an  alternate  set  of  equations  to  be  used  when 
singularity  conditions  are  approached  and,  equally  important,  the  means  to 
convert  to  this  set  during  flight  simulation. 

As  indicated  under  the  General  Equations  of  Motion,  special  attention  must 
be  given  to  the  equations  of  motion  when  ft  approaches  *  rr/2.  As  had  been  antic¬ 
ipated,  the  term  <,y  ..tan  ft  in  the  equations  of  morion  approaches  infinity  as  ft  ap¬ 
proaches  90  .  Furthermore,  0  becomes  indeterminate  when  this  singular  condi¬ 
tion  is  reached.  Unfortunately,  ft  90  occurs  whenever  the  missile  axis  be¬ 
comes  vertical,  the  occurrence  of  which  cannot  be  ignored.  To  avoid  this  condi¬ 
tion,  the  (H)  coordinate  system  will  be  replaced  by  a  new  fixed-plane  sys¬ 
tem  labelled  the  (VI  coordinates.  As  outlined  in  the  beginning  of  the  pro¬ 
cedure,  |y  is  constrained  to  tie  perpendicular  to  the  ife-  axis  (i.e.,  to  lie  in 

a  specified  fixed  vertical  plane)  while  iy  is  to  be  coincident  with  the  mis¬ 
sile  axis.  This  change  necessitates  redefining  the  Euler  angles,  and  deriv¬ 
ing  equations  oi  transformation  between  the  two  coordinate  systems. 

To  satisfy  the  conditions  imposed  upon  the  (V)  coordinate  system,  the 
rotations  indicated  below  are  made. 

1.  A  rotation  about  ij..  of  magnitude  4'. 

2.  A  rotation  about  j'  of  magnitude  so  that  «'  '  becomes  the 
missile  axis. 

To  avoid  ambiguity,  the  following  angle  restrictions  were  observed: 

0  V>0 

0  -ft  180 

In  matrix  form 


■U 


After  multiplying,  we  have 


*v 

C0 

ses'p 

-SGC'P 

‘E 

>v 

kv 

0 

S0 

-cosy 

S  Y 

C0CW 

i  e 
.  kE 

(88) 

Note  that,  as  required,  iv  has  no  component  along  iE. 


Rather  than  introduce  new  notation,  a  prime  over  a  given  "omega”  will 
denote  that  the  particular  angular  velocity  that  was  previously  referred  to 
in  (H)  coordinates  will  now  refer  to  the  (V)  coordinates. 


Continuing  in  a  manner  analogous  to  the  horizontal  coordinate  system, 
we  can  write 

«•  =  ’PT^  +  0“=’i’[C0l^  +  S0k^]  +  0j^  (89) 


or  iii  component  form 


and  for  later  purposes 


“kv  *  *» 


CO 


YV 


=  0 


(O 


ZV 


=  ¥S0 


CO 


CO 


ZV 


XV 


tan  0 


(90) 


(91) 


Note  that,  when  0  =  0°,  a  singular  condition  is  again  present;  however, 
investigating  the  geometry  of  the  rotations  for  V  and  0  will  show  that  this 
occurs  when  the  missile  axis  is  horizontal  rather  than  vertical,  a  situation 
for  which  the  previous  set  of  equations  is  applicable. 

The  angular  velocity  of  the  earth  in  the  new  coordinates  becomes 

Si  =  QiT*  =  [-ns0c^i“  +  ns^j”  -  ncoc^  (92) 
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Letting  ft'  be  the  angular  velocity  of  the  (M)  coordinates  relative  to 
(V)  coordinates,  we  can  write 


(wip  -  f!'  +  «'  d'  +  w'  +6 


(93) 


or,  in  component  form  (after  using  Equations  91  and  92), 


(<o'j 


ft' 


rxv  xvTtan© 


zv  -qsocv 


(w’ )  -  w'  +  QSV 

TYV  YV 


(94) 


KW ’  "kv  * ac0c'1' 


Differentiating  Equation  94  yields 


)  «  Q<  ^ZV  0  +  “zv  _  Q  [C0CTO  -  S0SVV] 

T  X  V  XV  sin2  0  tan  0 

(<4)yv  =  +  ncvv 

(o^zv  -  <^zv  -  Q  ts0c’y®  +  C0SV  vi 


(95) 


One  may  note  that  Equations  8  and  12,  the  basic  component  equations 
of  motion  (as  well  as  most  forces  and  moments),  were  derived  without 
first  specifying  the  orientation  of  the  (H)  coordinate  system;  hence,  they 
are  equally  valid  in  the  (V)  coordinates.  Therefore,  using  the  "omegas” 
and  the  "omega-dots”  derived  above,  along  with  the  new  transformation 
equations,  produces  an  alternate  set  of  equations  valid  for  vertical  orien 
tations  of  the  missile. 

Since  both  sets  of  equations  have  their  own  singularity  conditions,  it 
it  necessary  not  only  to  have  both  sets  of  equations  available  but  it 
should  also  be  possible  to  convert  from  one  system  to  the  other  as  the 
need  arises.  It  is  the  purpose  of  the  next  few  equations  to  establish  this 
conversion. 

By  definition  the  T  axes  of  the  (H)  and  (V)  coordinates  are  coincident; 
hence,  if  an  additional  rotation  is  made  about  the  iy  axis  so  th:<  jy 
rotated  lies  in  a  horizontal  plane  (as  determined  by  i^.  and  j^.  ),  then  the 
(V)  coordinate  system  will  become  coincident  with  the  (H)  coordinates. 
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Let  us  denote  this  additional  rotation  by  0,  whose  magnitude  is  yet  un¬ 
known.  In  matrix  form 

10  0  i  Y 

j~  =  0  CD  SD  j7  (96) 

k^  0  -S<1)  CD  kv 

Combining  the  inverses  of  Equations  88  and  96  eventually  produces 


C0 

S0SD 

S0CD 

S0SD 

CDCD  -  C0SDSD 

-C'VSD  - 

-S0CD 

S^CD  f C0CDSD 

-S*PSD  4 

Comparing  this  with  Equation  17  yields  explicit  expressions  for  ob¬ 
taining  V,  0,  and  D  in  terms  o(  6  and  0 ,  as  follows: 

Cos  0  cos  0  cos  0  (98a 

Sin  D  -  cos  (9  sin  jj .  cos  ip  sin  6  (98b 

sin  0  sin  0 


Sin  D  -  ~siri  V 
sin  0 


cos  D  ^  sin  6  cos 
sin  0 


Again  no  ambiguity  results  from  the  restrictions  made  earlier  on  0  (i.e., 
the  quadrants  of  0,  V,  and  D  are  uniquely  determined). 

Knowing  S',  0,  and  D  and,  of  course,  0,  0,  0,  and  (),  we  are  now  in  a 
position  to  determine  the  angular  rates  cf  change,  H*  and  0. 

This  is  accomplished  by  differentiating  Equations  98a  and  98b,  producing 

<3  SflC0< b  +  CAS00 

S0 


^  sw|  -sos^o  *  oc00 1  -  a)S0c00 

CV  sin'  0 


if  ¥  happens  to  be  either  90°  or  270°,  one  can  differentiate  the  second 
equation  (98b),  obtaining  instead  of  Equation  100 


^  =  -s0c<9<9  +  sflcee 

SW  sin1© 


(10  1) 


to  eliminate  the  singularity. 

Once  these  derivatives  are  known  it  follows,  from  Equation  90,  that 


o>2V  =  V  sin  0 


6)'  =0 
YV 


(102) 


To  obtain  Oxv,  one  need  only  equate  the  two  expressions  for  (wT)x, 
obtaining 


^  .  -  a  S 0C¥  =  flxH  -  <uZH  tan  Q  -  Q  SS  ( 103) 
tan  0 


or,  since  S0CW  ^  S0, 


fiXV  "  nXH~wZHtan^' 


to 


ZV 


tan 


0 


(104) 


To  complete  the  transformation,  Equation  96  again  comes  into  use.  We 
have  for  the  invariant  vectors  (force,  moment,  and  velocity) 


F' 

XV 

-  F 
rXH 

LXV 

=  lxh 

F' 

YV 

=  F  C4>  -  F 
YH  rZH 

S«t> 

4v 

=  lyhc<I>  “  lzhs4> 

FL 

=  F  S<t>  +  F„, 

C<1> 

L’ 

=  L„,  S$  +  L  ,UC<1> 

ZV 

YH  ZH 

ZV 

YH  ZH 

V' 

-  V 

XV 

XH 

V' 

-  V  C4> 

-  V 

S<U 

YV 

YH 

ZH 

V' 

=  V  S<t> 

4  V 

CO) 

ZV 

YH 

ZH 

(105) 
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Other  equations  required  for  using  the  fixed  vertical  plane  are 


a;,  -  -Sac©  f  Casesv 

A|,  S/}CaCfl  .  S£S«S0S>f'  -  CftSQC'V 


(108) 


This  result  and  the  other  two  components  are  tabulated  in  the  Summary  of 
Equations  (p  51).  Needless  to  say,  (V  )  is  determined  from  Equation  72 

written  in  terms  of  the  (V)  coordinates.  Wind  for  the  flat  earth  case  is 
analogous  to  that  of  the  (H)  coordinates.  The  resulting  equations  are  also 
tabulated  in  the  Summary  of  Equations  section. 


Before  leaving  this  section,  three  derivations  are  yet  required.  These 
involve  expressing  the  thrust  and  thrust  moments,  specifying  new  initial 
conditions  in  the  (V)  coordinates,  and  converting  from  (V)  back  to  (H). 
We  begin  with  the  thrust  equations. 


Since  $  is  defined  as  the  angle  between  the  j  and  j  axis,  Figure  8 

V  H 

(p  26)  may  be  modified  as  follows: 


4c 


from  which  it  is  clear  chat  the  thrust  components  expressed  in  the  new  co¬ 
ordinates  are  given  by  Equation  58  with  S  replaced  by  S' .  Thus, 

Txv  =  T  cos  ST 

TyV  =  T  sin  8tCos(§a+  S')  ( 109) 

Tzv  =  T  sin  §T  sin  (5A  +  S') 

At  changeover  (i.e. ,  when  converting  from  (H)  to  (V)  coordinates),  angle 
S'  is  defined  by 


S'  =  <t>  +  S 


(110) 


after  which  the  following  is  used 

S’  =  f*  O'dt  +  S'  I 

t  changeover  It  changeover 


Use  of  these  equations  implies  that  there  is  no  reorientation  of  the  (M) 
coordinates;  hence,  all  quantities  previously  referred  to  these  coordinates 
(such  as  SA,  Sp  rXM’  rYM>  rZM^  retnain  constant  when  changing  from  (H) 

to  (V)  coordinates. 

One  similarly  replaces  S  by  S'  in  Equation  62  to  obtain  the  thrust  mo¬ 
ments  in  the  (V)  coordinates. 

This  completes  the  transformation  equations  pertinent  to  the  thrust  terms. 
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Initial  conditions  must  be  provided  in  the  (V)  coordinates,  because  of 
the  possibility  of  "vertical  launch."  Actually,  much  that  has  already  been 
derived  in  the  section  on  initial  conditions  can  be  used  here.  The  only 
real  change  is  that  angle  L  has  to  be  redefined  relative  to  the  (V)  coordi¬ 
nates.  Angles  A,  B,  G,  and  H  are  defined  as  before.  Let  L'  denote  this 
new  angle.  It  is  determined  by  noting  that  jv  is  to  have  no  component 

along  iE  ,  in  place  of  kp.  This  is  effected  by  equating  Equation  80  (with 
H’s  replaced  by  V’s)  with  Equation  88  (instead  of  with  Equation  16  as 
was  done  before).  Mathematically,  we  have 

Cos  0(0)  -  bu 

Sf  C¥  -JjiL  (112) 

s?T  S0 


We  have  for  L' 


CL* 

SL' 


b„  CV 

bj  i  b,j  -  bj,  bj } 

b3,CW 

b:  i  bjj  —  b,,  bj , 


(113) 


where  blt,  ...,  b,,  are  defined  as  before. 

For  the  flat  earth  case,  again  setting  A  0°,  B  270°,  we  have 


Cos  0(0)  SHCG 

Cos  V(0)  -Q1  (114) 

S0 

Sin  M'(0)  -  -§G 

S0 


Finally,  it  is  necessary  to  determine  expressions  for  converting  from 
(V)  to  (H)  coordinates.  To  obtain  these  expressions,  one  can  proceed  in  a 
manner  similar  to  previous  derivations. 

An  angle  f  is  defined  such  that  when  a  rotation  of  this  magnitude  is 
performed  about  iv,  the  jv  axis  will  then  lie  in  the  horizontal  plane  (i.e., 
parallel  to  the  equatorial  plane).  This  rotation  expressed  in  matrix  form 
will  yield  Equation  96,  with  $  replaced  by  <t>'  . 
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Postmultiply ing  the  inverse  of  Equation  88  by  the  inverse  of  Equa 
cion  96  produces 


1 

- 

~ 

‘K 

c:0  s0s<t>‘  s(-k  <h’ 

*H 

•e 

30S4»  efeth’  -  CfisVsfl*'  -CH'Sth’  -  COSl'C^' 

ill 

k  E 

-SfiC'f'  SH'C  <ti'  ,  C0CYS<D'  -S4>'  SV  +  C0C<fc'(  V 

.  kH. 

(115) 


Comparing  Equation  115  with  Equation  17  below 


coco 

-Sy) 

soco 

easy) 

co 

soso 

-so 

0 

CO 

results  in  the  following  transformations: 


SO  S0C¥ 


CO 

C0 

sc- 

S0S¥ 

CO 

co 

s<h' 

-SO 

c<t>' 

soco 

S0 

S0 

(17) 


(116) 


No  ambiguity  results,  since  -1  _o  •. 


To  obtain  the  angular  velocities,  one  differentiates  Equation  116, 
obtaining 


-)  CHCVH--SHSW 

CO 

CO  cos'  0 


If  0  is  90  or  770  ’,  one  can  use 


u>s00  -  c.Qsoo 

Si,'/  cos'1  0 


(117) 


(118) 


(119) 


Equation  104  completes  the  transformation  of  the  angular  velocities. 
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Solving  for  QxH^eac^s  to 


fiXH  Oxv  +  w7H  tan  ^  t  — — —  ( 120) 

tan  0 

where  both  oy  M  and  are  defined  by  Equation  23- 

When  Equation  96  is  used,  the  forces,  moments,  and  velocities  can  be  ex 
pressed  in  (H)  coordinates,  given  their  representations  in  (V)  coordinates. 
These  auxiliary  equations  are  tabulated  in  the  following  section  of  this 
report. 

This,  therefore,  completes  the  transformation  from  the  (V)  to  the  (H) 
system.  A  complete  tabulation  is  presented  in  the  following  section. 

Summary  of  Equations 

The  equations  of  motion  as  derived  in  the  preceding  five  sections  are 
summarized  below.  In  the  writing  of  the  equations,  the  primes  have  been 
dropped  when  working  in  the  (V)  system.  No  ambiguity  should  result. 


SO 


Force  Equotions  (H-Coordinotes) 
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tan  6  -  8S0I  ly  wy(l  -  rL»Y(J  Ix  “  <'>z„  tan  ^  “  8S011 


where 
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II.  For  The  Flot  Earth  Cose 
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where  if 
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Equations  of  Transformation 
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SHCBSA  -  CHSGCA  +  CHCGSBSA 


CHSB  -  SHCGCB 
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=  0  if  trajectory  starts  in  H-Coordinates,  T  =  Tchangeover  otherwise. 


In  all  equations 
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0  JOD  ^Z(|)  _  Q  ubj  HZ (7)  _  HX; 


Fore*  Equations  (V-Coordinates) 
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— pd2  (Vj.)2  kjg  V  '  (V  )yy  (r^-.  -  A  jj)  —  pd  ( \  )  kp  V  (V  )^y  (Xq  A^)  pd4  (Vf)  k 
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ce  sy  +  rze  cec4» 


||.  For  The  Flat  Earth  Cose 
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=  (Vw)yw  a'ji  +  <vw>zw  a' 
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(vr)xv  =  [vXE  -<vw)xFlce  +  [vYE  -(VW>YE1S©SV-[VZE-(VW)ZE]S0CSI 


=  [VYE  -(V^glCf+hfeE  -(VW)ZE1SV 

=  lVXE  -(Vw)xE]S0-[VYE  -(Vw)Yl£]C0SV+[VZE  -(Vw)ZE]C0CS» 
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0  if  trajectory  starts  in  V-Coordinates,  T  =  ^changeover 


Equations  of  Tronsfof motion  (V)  -*  (H) 
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In  all  equations 


RESULTS  AND  DISCUSSION 


The  final  equations  resulting  from  the  derivation  are  tabulated  in  the 
preceding  section.  There  are,  however,  many  aspects  of  the  equations  that 
have  not  been  discussed.  It  is  the  purpose  of  this  section  to  consider 
briefly  some  of  these  aspects  and  limitations  and  to  mention  what  is  in¬ 
volved  in  incorporating  certain  refinements  in  the  equations  as  presently 
derived. 

One  might  begin  by  noting  that  a  typical  rocket  trajectory  consists  of 
several  flight  phases.  Explicitly,  these  phases,  typical  for  a  two-stage 
rocket,  may  be  tabulated  as  follows: 

1.  Acceleration  of  booster  and  main  stage 

2.  Coasting  of  booster  and  main  stage 

3.  Separation  of  booster  and  main  stage 

4.  Coasting  of  main  stage 

5.  Acceleration  of  main  stage 

6.  Free  flight  of  main  stage. 

These  phases  are  in  direct  correspondence  with  those  used  in  the  six- 
degree-of-freedom  digital  computer  trajectory  program  at  Picatinny  Arsenal. 
The  user  is  thus  able  to  run  any  or  all  of  these  phases.  Phase  3  is  effected 
by  imparting  a  short  thrust  to  the  main  stage,  simulating  the  actual  booster 
separation. 

Since  the  aerodynamic  coefficients,  ccnte«s  of  gravity,  pressure,  etc.  are 
not  necessarily  identical  for  all  phases,  several  sets  of  such  data  are 
often  required  to  simulate  a  complete  trajectory.  The  equations  in  the  pre¬ 
ceding  section  of  thL  report  a.e  used  tor  ail  phases. 

To  numerically  solve  the  differential  equations  of  motion,  a  modification 
of  the  third  order  Runge-Kutta  technique  is  used  (Runge-Kur  >Gill).  This 
numerical  scheme  appears  to  be  adequate  for  both  low  -  and  high-spin  pro¬ 
jectiles,  although  inclusion  of  thrust  malalignments  remains  to  bt  investi¬ 
gated.  As  is  indicated  in  the  introduction,  however,  i:  appeared  that  th;s 
numerical  scheme  was  not  satisfactory  when  high  angular  'ates  of  change 
were  implicit  in  the  equations  themselves.  This  condition  pertains  to  high- 
spin  projectiles  for  which  the  fixed-plane  coordinate  system  is  rigidly 
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attached  to  the  missile  (as  was  originally  done).  This  is  an  advantage  of 
the  present  coordinate  system. 

Other  aspects  of  these  equations,  such  as  the  use  of  thrust  modification 
factors,  the  atmospheric  model  (ARDC  Atmosphere  of  1959).  etc.,  will  be 
discussed  in  a  forthcoming  report. 

Some  of  the  refinements  of  the  equations  are  discussed  below. 

1.  Introduction  of  an  ellipsoidal  earth  model.  This  modification  has 
far-reaching  effects  on  the  equations.  For  example,  the  direction  and  mag¬ 
nitude  of  the  gravitational  attraction  (consistent  with  the  earth  model) 
would  vary  with  missile  position.  Additional  coordinates  would  have  to  be 
specified  on  the  surface  of  the  ellipsoid  to  properly  introduce  wind  data. 

In  addition,  an  iterative  scheme  is  necessary  to  compute  the  altitude  of  the 
missile  if  it  is  naturally  defined  as  the  shortest  distance  from  the  missile 
CG  to  the  earth  model.  Initial  conditions  would  also  be  modified  accordingly. 
For  increased  ranges  and  accuracy  requirements  of  trajectory  simulation, 
the  ellipsoidal  earth  refinement  may  well  become  necessary. 

2.  Treatment  of  asymmetrical  missiles.  Asymmetrical  missiles  can 
arise  from  two  sources.  The  more  severe  situation  is  present  if  the  ex¬ 
ternal  missile  configuration  does  not  possess  rotational  symmetry.  For 
this  case,  the  transverse  moments  of  inertia  could  vary  about  the  fixed- 
plane  coordinates  as  the  missile  rotates.  Also,  specific  data  relating  the 
resultant  point  of  application  of  the  aerodynamic  forces  would  have  to  be 
provided  for  accurate  results.  A  second  asymmetrical  condition  can  occur 
for  externally  symmetrical  missiles  when  the  missile  CG  is  offset  by  a 
prescribed  amount  from  the  longitudinal  axis  of  the  missile.  This  obviously 
modifies  the  lever  arms  in  the  moment  equations  and  may  also  introduce 
time-dependent  moments  and  products  of  inertia. 

3.  Other  refinements  could  include  guidance  and  launcher  effects,  both 
of  which  can  play  prominent  roles  in  trajectory  analysis. 

In  addition  to  the  limitations  listed  earlier,  two  assumptions  were  im¬ 
plicit  in  the  derivation.  The  first  is  that  the  motion  of  the  earth  about  the 
sun  was  neglected;  for  most  sub-orbital  trajectories  this  phenomenon  can 
be  ignored.  The  second  assumption  is  that  the  rate  of  change  of  inertias 
was  neglected  (i.e.,  i  terms  during  the  thrusting  periods).  For  missiles 
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possessing  excessively  high  burning  rates  with  large  angular  velocity  com¬ 
ponents,  the  magnitude  of  these  terms  should  be  investigated. 

A  final  word  about  the  equations  is  that  there  is  no  estimate  of  the  dis¬ 
persion  of  the  missile.  This  requires  the  computation  of  several  trajectories, 
each  for  a  slightly  different  initial  condition,  with  appropriate  statistical 
combinations  of  the  various  ranges  and  deflections  from  a  fixed  standard. 
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APPENDIX  A 


Rotating  coordinate  systems 

It  is  the  purpose  of  this  Appendix  to  express  the  rate  of  change  of  a 
vector  relative  to  fixed  coordinates,  in  terms  of  rotating  coordinates.  Briefly, 
this  analysis  answers  the  following  question:  "Given  a  particle  whose 
motion  is  known,  how  would  the  particle’s  motion  appear  to  someone 
situated  on  a  coordinate  system  that  is  itself  in  rotation?” 

Obviously,  one  must  know  the  motion  of  the  rotating  coordinates  relative 
to  the  fixed  coordinates.  To  be  explicit,  let  (I)  denote  fixed  coordinates 
and  (H)  a  rotating  coordinate  system,  whose  origins  are  both  coincident. 
Further,  let  <ibe  the  angular  velocity  of  (H)  relative  to  (I).  This  motion 
might  be  as  depicted  in  Figure  A-l,  where  the  three  dotted  discs  represent 
the  paths  of  each  axis  of  the  (H)  system  for  the  given  constant  u. 


Fig  A-l 
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Let  Rbe  an  arbitrary  vector.  Clearly,  R  has  a  representation  in  the  (H) 
system  as  follows: 


R  =  RXH  ‘H  +  RYH  >H  +  RZH  kH 


(A-l) 


To  determine  the  rate  of  change  of  ft  relative  to  the  (I)  coordinates  in  terms 
of  the  (H)  coordinates  one  must  account  not  only  for  the  changing  magnitude 
of  ft,  but  also  for  the  variation  of  the  unit  vectors  iH,  jH,  kH  (which  help 

represent  ft)  relative  to  (I).  Mathematically: 


UiA  ,  , 

—  =  rxh  ‘h  +  ryh  >H  +  RZH  kH  +  rxh  ~r 
dt  dt 


dt 


dt 

(A-2) 


The  first  three  terms  of  Equation  A-2  define  simply  how  ft  itself  is  changing 
independent  of  any  moving  coordinates.  The  latter  three  terms  describe  the 
motion  of  (H)  with  respect  to  (I).  It  is  noteworthy  that  the  I  subscript  de¬ 
notes  differentiation  with  respect  to  the  (I)  coordinates.  To  elaborate  fur¬ 
ther,  if  I  were  replaced  by  H  then  Equation  A-2  would  read 


dt 


RXH  ‘h 


ryh  ’h 


RZHkH 


(A-3) 


since,  for  purposes  of  the  differentiation,  iH,  jH,  kH  would  be  fixed. 

•  ▼  d  lc 

The  task  remains  of  obtaining  expressions  for  .  .  .  ,  ““jp*  To  do 

this,  let  us  "isolate”  the  i  vector  as  presented  in  Figure  A-l,  and  define 
some  new  quantities  as  given  in  Figure  A-2.  Further,  let 
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r.  . 


iH  (t  +  At)  denotes  the  position  of  iH  at  time  t  +  At.  Clearly,  by  the  defi¬ 
nition  of  a  derivative, 


iH  (t  +  At)  -  iH  (t) 


Liin  .  . 

dt  At-»o  At 


and  from  the  geometry  of  Figure  A-2 

iH  (t  +  At)  -  iH  (t)  =  A  iH  ~  iH  sin  0  [«At] 

we  may  write 


(A-4) 


(A-5) 


dI  ‘H 

i^  sin  0  « At 

dt 

At 

and  in  the  limit 


jljjj 

dt 


=  Lim 
At-0 


iH  sin  0  oAt 

At 


iH  sin  0 M 


(A-6) 


Equation  A-6,  however,  is  very  reminiscent  of  the  cross-product,  namely 

(A-7) 


dirH 


dt 


=  MX  1 


lh 


and,  if  one  examines  the  directions,  one  can  indeed  state 
dl*H  -  T 


dt 


*mx  iH 


(A-8) 


Since  iH  is  typically  representative  of  each  coordinate  axis,  we  also  have 


*H  -  r 

—  a  M  X 


di  ku  ■*  ~ 

JL”  =  MX  L 

dt  ^ 


(A-9) 


71 


It  is  natural  to  substitute  Equations  A-8  and  A-9  into  Equation  A-2.  This 
produces 


djR  •  —  •  —  *  —  -»  — 

~7~  =  rxh  ‘h  +  ryh  *h  +  rzh  ^  +  rxh  (MX 


+  Ryu  (wx  iH)+  Rzh  (“x  kH^ 


which  one  may  write  as 


^— =  ^4—  +  wx  rxh  ‘h  +  “  x  ryh  Jh  +  “*  Rzh  kH 

dt  dt 

or  finally 

d.R  d„R  -  -* 

— ! —  =  — —  +  u  x  R 

dt  dt 

which,  in  the  text  proper,  is  the  basis  for  Equation  4. 


(A- 10) 


(A- 11) 


( A  - 12) 
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APPENDIX  B 

Matrix  representation  of  rotations 

Assume  the  geometry  of  Figure  B-l,  where  a  rotation  ofmagnitude  0  has 
been  performed  about  the  k£  axis,  producing  new  vectors  i'  and  j’ . 


Fig  B-l 


Clearly  the  projection  of  i'  on  the  iE  axis  is  i'  cos  0  =  cos  0 

(all  vectors  shown  are  of  unit  magnitude).  Continuing  in  like  manner  for 

all  possible  combinations,  one  may  form  the  following  table: 

‘i  >1 

i*  cos  0  sin  0  0 

j*  -sin  0  cos  0  0 

k'  0  0  1 

Here,  for  example,  -sin  0  is  to  be  interpreted  as  being  the  projection  of  j' 
on  the  ie  axis. 

If  we  now  choose  to  represent  i'  in  terms  of  its  components  along  i  ,  j  , 
and  k£,  we  write 

i*  *  cos  0  iE  ♦  sin  0  jg  ■»  0  kg  (B-  1) 
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Though  matrix  multiplication,  we  can  write  not  only  i',  but  also  j*  and 
k' ,  not  separately  as  in  Equation  B-l,  but  together  in  one  matrix  similar 
to  the  given  table.  Thus; 


i' 

cos  0  sin  0  0 

’  »E 

i' 

= 

-  sin  0  cos  0  0 

*E 

k' 

0  0  1 

.y 

or  Equation  13  of  the  text. 

It  might  be  mentioned  that  this  development  is  not  intended  to  be  rigorous; 
however,  it  should  enable  one  to  obtain  these  rotation  matrices  without 
difficulty. 

Finally,  let  us  complicate  the  above  and  now  assume  that  i'  and  j*  are 
rotating  about  kg  with  some  angular  velocity  u,.  Clearly  u,  is  directed 
along  the  kg.  axis,  and  further,  in  time  At  we  have  the  relation 


i«,i  ■ 


(B-3) 


or,  in  the  limit,  as  At  approaches  zero 


!w,l 


dt 


(B-4) 


Finally,  combining  both  magnitude  and  direction,  we  can  write 

«•  (B*5) 

Ve  can  do  likewise  for  a  rotation  about  j' .  obtaining,  for  example. 


«*,  « i'  (B-6) 

where  9  is  an  angle  defined  analogously  in  the  i*  -  k*  plane  as  0  was  de¬ 
fined  in  the  ie  -  jg  plane. 

Since  angular  velocities  may  be  added,  we  can  combine  Equations  B-3 
and  B-6  and  write 

ii*'  • 

tit  in |  ♦  hi i  0  kg  •  ti  j'  (B-7) 

which,  in  the  text,  is  Equation  22. 
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